Are the Tails of Percolation Thresholds Gaussians 
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Abstract: The probability distribution of percolation thresholds in finite 
lattices were first believed to follow a normal Gaussian behaviour. With 
increasing computer power and more efficient simulational techniques, this 
belief turned to a stretched exponential behaviour, instead. Here, based on 
a further improvement of Monte Carlo data, we show evidences that this 
question is not yet answered at all. 

In reference |T], the percolation on a A^-site square lattice is treated with 
high numerical accuracy. Indeed, the best known estimate for the critical 
threshold, pc = 0.59274621(13), comes from this work. In order to study this 
kind of problems, the authors follow a very fruitful Monte Carlo approach 
which allows one to obtain continuous functions of p, the concentration of 
occupied sites, namely the canonical-like average 

m = E p"(i-p)^-"i?„ , (1) 

n 

of some quantity R. Here, i?„ is a uniform average over all configurations 
with just n occupied sites, i.e. a microcanonical-like average. are simple 
binomial factors. By filling the initially empty lattice, site by site at random, 
and repeating this process many times, one is able to get the discrete set of 
microcanonical averages i?„ accumulated into an ra-histogram, over the entire 
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Figure 1: Spanning probability function for a fixed lattice size L (solid line). 
For larger and larger sizes (dotted lines), this function approaches a step. By 
fixing some value r at the vertical axis, one can find a sequence of values piif) 
at the horizontal axis approaching the critical threshold Pc, for increasing 
lattice sizes. 

range, n = 0, 1, 2 . . . A^. From this set of numbers, the determination of the 
continuous p-function R{p) is straightforward. 

In particular the authors of ^ fix attention at the horizontal wrapping 
probability Rl{.p) around a. L x L torus, i.e. a square lattice with periodic 
boundary conditions on both directions. In the thermodynamic limit, this 
function approaches a step: -Roo(p) = below the critical threshold Pc, and 
Rooip) = 1 above pc- For finite sizes, Rl{p) presents a sigmoid aspect similar 
to figure 1. A good approach to pc is to choose some fixed value r, and solve 
the equation Rl{p) = r, getting the root p shown at the horizontal axis. 
Here, one can appreciate the advantage of knowing Rl{p) as a continuous 
function of p. Keeping the same value r and repeating this task for a series of 
increasing lattice sizes (dotted lines), one gets a series of roots ^^^(r), PljI'")' 
PL-ji'r) ■ ■ ■ which converges to the desired threshold pc in the thermodynamic 
limit. 

The above reasoning is valid no matter which is the fixed value for r 
one chooses. However, for the very particular choice r* = 0.521058290, 
a universal probability exactly known through conformal invariance argu- 
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ments j2], the convergence becomes fast, i.e. the root p{L) differs from pc as 
j^-2-i/u _ /^-2.75^ where u = 4/3 is the correlation length critical exponent. 
The above quoted accurate value for pc was obtained in this way. For details, 
see and references therein. 

Reference j2] proposes the mathematical form 



for estimators pi obtained from quantities like R{p). The option for the 
wrapping probabilities around the torus and the convenient choice of Pinson's 
number r = r* lead to vanishing values for the two first terms Ao(r*) = 
Ai{r*) = 0, a lucky coincidence which accelerates very much the convergence. 

The quantity Rl{p) is obtained, as quoted before, by filling up the initially 
empty lattice site by site, at random. Clusters of neighbouring occupied sites 
grow. As soon as the horizontal wrapping along the torus is set, one books 
the corresponding value of n, the number of occupied sites so far, and stops 
the process. For that particular sample, the wrapping probability is a 
step function, i.e. Rl = below n and Rl = 1 above. The same routine 
is repeated many times, in order to have a probability distribution for n. 
The various step functions are then superimposed to get the microcanonical 
averages Rn in equation (1), stored in an n-histogram. Finally, the continuous 
canonical average Rl{p) can be calculated for any value of p. 

Each process of filling-up the lattice (one sample) yields a single value 
n for the statistics, i.e. just one more entry on the n-histogram. In jB], 
we decided to improve this point, by changing the definition from wrapping 
to spanning probability, figure 1. We fix two parallel horizontal lines on 
the L X L torus, separated by a distance of L/2, for instance lines i = 1 
and i = 1 + L/2. The measured quantity is now the probability of having 
these two lines connected by the same cluster of neighbouring occupied sites, 
instead of the wrapping probability along the whole torus. The advantage 
is that we can measure the same thing for lines i = 2 and i = 2 + L/2, also 
for lines i = 3 and i = 3 + L/2, and so on. Moreover, also vertical parallel 
lines can be included into this counting. At the end, from a single sample 
we store L new entries into our n-histogram, instead of just one more entry. 
Note that this advantage even increases for larger and larger lattices. 

Within the same computational effort, our approach allows the test of 
larger lattices. Because of that, we were able to confirm the validity of 
equation (2) with high precision, by sampling 27 different lattice sizes from 
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L = 18 up to L = 1594, a 8-thousand factor in number of sites. On the 
other hand, our definition does not allow the chance of vanishing both Ao(r) 
and Ai{r) at once. We can have only Ao{f) = 0, for a particular universal 
probability f = 0.984786(11) numerically determined within the same work 
jS]. Independently, Cardy |3] tried to determine it by conformal invariance 
arguments, however, in looking for configurations which link two parallel 
lines, he was forced to disregard configurations which also wrap along the 
other direction. As a result of using larger lattices but a slower convergence 
rate of L"^"^/^ = L'^ '^^, we get the same figure Pc = 0.59274621(33) as 
in jT], but within a 3 times larger error bar. 
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Figure 2: Test of equations (3) or (4), for tails on the left of figure 1. The 
five continuous curves correspond to L = 1594, 1354, 1126, 958 and 802, 
from left to right. In each case we sampled 4 million lattice-filling processes, 
which corresponds to 6, 5, 5, 4, and 3 xlO^ entries in each n-histogram, 
respectively. The statistics is improved by a factor over than 1000, compared 
with pp for equivalent lattice sizes. The dashed lines show the alternative 
slopes 2 (right) or 4/3 (left). In the authors' opinion, no definitive conclusion 
is possible. 

The non-Gaussian behaviour of the finite-lattice-threshold distribution 
near the infinite-lattice critical point is already stablished Here, we profit 
from the same simulational data in order to investigate the distribution tails, 
far from the critical point. Which is the mathematical form of the tail 
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Figure 3: Alternative test of equations (3) or (4), for the same lattice sizes 
L = 1594, 1354, 1126, 958 and 802, from left to right. 

observed in figure 1, below the root p? One possible answer is a simple 
Gaussian form [HI Ej 

RL{p)-^exp[-Kip-p,f] . (3) 
Another alternative is a stretched exponential jHl El QD] 

RLip) = eM-CiPc-pr] , (4) 

where the strict equality (for large L and far from p^) is a consequence of the 
periodic boundary condition [1 , which holds for our data. Profiting from this 
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strict equality, one can test equation (4) by constructing a plot of ln[— ln(i?)] 
against \n{pc — p). This was done in P, and we repeat the same for our 
data, in figure 2. Note that our range for ln(pc —p) (up to —3.4 for L ~ 10^) 
is larger than in reference PP (up to —4.6 for the same size). This means 
that we are testing more deeply the distribution tails, thanks to our trick of 
sampling L new entries for each run. Even so, the conclusion in favour of 
either equation (3) or (4) is by no means obvious. Note a further difficulty 
in what concerns equation (3), because the leading multiplicative constant 
in front of the exponential is not necessarily 1. 

Another, perhaps better way to address the same question is by plotting 
ln(i?) twice, against (p—pc)'^ and {pc—pY^^- Figure 3 shows the result for our 
data. Note that our range for ln(i?) (down to —16) doubles the one presented 
in [T]. The would-be Gaussian case (up) presents clear positive curvatures, 
whereas the would-be stretched exponential (down) presents negative curva- 
tures although not so pronounced. The exponent 4/3 seems to fit better, but 
one cannot extract a clear conclusion from these data. 

Still more undefined is the situation of the other tails on the right of figure 
1, above Pc- In this case (not shown), our accuracy limit for ln(l — R) (down 
to —16) is reached much closer to Pc than the case shown in figures 2 and 3, 
below Pc- 

Concluding, we present new Monte Carlo data concerning the distribution 
probability of percolation thresholds on a finite square lattice. We address 
the question of the mathematical form of the distribution tails, equation (3) 
against (4). Even considering that our statistics is over than 1000 times 
larger than previous works, no definitive conclusion can be extracted from 
our data, in what concerns the asymptotic tail exponent. 
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